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Abstract 

We study the behavior of metric perturbations in a recently proposed model of 
phantom dark energy with tachyonic instability at long wavelengths. We find that 
metric perturbations exponentially grow in time, starting from very small values deter- 
mined by vacuum fluctuations, and may become sizeable at late times. This property 
may be of interest for phenomenology. 

1 Introduction 

Among models of dark energy (for reviews see, e.g., Refs. [TJ, El El HI El El E] ) ) perhaps the 
most exotic ones are those exhibiting phantom equation of state, p = wp with w < —1. 
Observational data strongly constrain this possibility [SJ El [101 [HI [121 03] , but by no means 
rule it out, especially if the equation of state parameter w varies in time (which is the 
case in many phantom models). For phantom matter, weak energy condition is violated, 
signaling instabilities and possibly other pathologies. Indeed, the simplest phantom models 
like those involving negative energy scalar field [H] have ghosts in the spectrum at arbitrarily 
high spatial momenta, and hence are plagued by catastrophic vacuum instability unless one 
introduces Lorentz- violating cutoff into the theory [T51 [TB] . 

This situation is not generic, however. First, phantom behavior may be mimicked in 
scalar-tensor and f(R) gravities [TTICEEICEH] which need not have instabilities at all. Second, 
in theories with Lorentz-violating phantom fields, the violation of weak energy condition at 
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cosmological distances does not imply dangerous instabilities at much shorter length scales. 
Indeed, phenomenologically acceptable models have been proposed [201 ED [22J [23] that have 
precisely this property. In particular, the model of Refs. [221 123] has tachyon in the spectrum 
of small perturbations about homogeneous Lorentz-violating phantom background, but the 
instabilities occur at low spatial momenta only, so that their time scale may be roughly 
comparable to (though somewhat smaller than) the present age of the Universe. 

The latter property may well be generic to a subclass of phantom models, in which 
phantom instability is of tachyonic nature. So, it is of interest to understand, using the 
model of Refs. [22;, 23J as a prototype, what observable consequences the long-wavelength 
tachyonic instabilities may have. As a step in this direction, we calculate in this paper the 
metric perturbations generated by the tachyon modes in this model. We find no surprise: 
like the tachyon modes themselves, metric perturbations exponentially grow in time, starting 
from very small values determined by the quantum physics of vacuum fluctuations. This may 
give rise to potentially observable effects, as the Universe rapidly becomes less homogeneous 
at late times. 

2 The model 

The model for dark energy we study in this paper is the Einstein gravity interacting with a 
vector field and a scalar field 0. Its construction combines two ideas. One [23] is that 
in the absence of gauge invariance, a theory of a vector field may nevertheless have healthy 
spectrum, provided that the background vector field is non-trivial. The other [22] has to do 
with the fact that with vector field(s), one can construct generally covariant one-derivative 
terms in the Lagrangian, that dominate at low momenta and become negligible at high 
momenta. Hence, the Lagrangian for the vector and scalar fields is [22] 

L= -l a (£)I^iV^ + ^^ , 

where £ = B^B^/A 2 , A is an UV cut-off scale, dimensionless parameters a and /3 are functions 
of £, and e is a positive constant of dimension of mass. Upon introducing the notation 
B 2 = B^B^, this Lagrangian can be rewritten as follows, 

L = -p fA D^ + Z+l W*)^*) + + *JB> + £b> - =V , (!) 

where = x5-/3(0 — To simplify formulas, in what follows we assume that = 

const, 7(£) = const. The range of parameters in which the model exhibits phantom behavior 
is 

M < m < e , (2) 
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and a ~ 7. Spatially flat homogeneous background in this model is described by 

g^v = a 2 (v)v^ , B o = a(rj)X(r]) , = ip{rj) , 

where rj is conformal time. Importantly, temporal component of the vector field does not 
vanish, which ensures that the expression ([1]) makes sense, and that the spectrum of per- 
turbations about this background is healthy (22J, [23] at high spatial momenta (no ghosts, 
tachyons or superluminally propagating modes), provided that 

a > 7 > . 

As we will see shortly, one of the scalar modes of perturbations becomes a tachyon at 
sufficiently low momenta. 

Cosmology in this model, with the usual and dark matter added, is fairly interesting [23J. 
In a wide range of initial data for X and (p, these fields stay constant at early, matter 
dominated stage, and the corresponding equation of state is w = — 1. Then they start 
evolving, the dark energy equation of state first being normal, w > — 1, and later phantom, 
w < — 1. The evolution often occurs in the slow roll regime, so that \w + 1| is small in the 
entire course of the cosmological expansion. Finally, the system approaches the de Sitter 
attractor, which to the leading order in m/e, M/e is characterized by the following values of 
the Hubble parameter and fields, 

M mM PL MM PL 

Ha = ' Xa = 7tt^ ' i f >A = n — • 

V3a ylZiie ■s/Anam 

The interesting regime in which dark energy contributes substantially to the total energy 
density, occurs at 

H ~ H A , X ~ X A , Lp ~ if A . (3) 
In what follows, we will be interested precisely in this range of the cosmological variables. 



3 Tachyonic perturbations 

The linearized perturbations of dark energy fields break up into vector and scalar parts. The 
vector sector is healthy at all spatial momenta, whereas the scalar sector has tachyons at 
relatively low momenta. Hence, we concentrate on the scalar sector. In Minkowski space- 
time, the momenta at which one of the modes is tachyonic are of order P ~ e, and the 
corresponding "frequencies" are of the same order p2, 23]. In the expanding Universe, 
a mode of a given conformal momentum p becomes tachyonic as the physical momentum 
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redshifts down to P = p/a ~ e. Our purpose is to calculate the metric perturbations 
generated by the tachyonic modest- 

With scalar perturbations included, the fields and metric in conformal Newtonian gauge 
are given by 

V177 X) 

B (t],x) = a{rj)X{rj) + &o(»7,x) , B^-q, x) = ^(77, x) , 0(r?,x) = ^(77) + 



0(77) 

ds 2 = a 2 (r]) [(1 + 2<&)dif - (1 - 2$)dx 2 ] 

where 6j is a longitudinal vector, and we have made use of the fact that the linearized 
energy-momentum tensor of the fields B^ and <fi has zero anisotropic stress. 

The complete expressions for the quadratic action and linear equations for perturbations 
in the slow roll regime (i.e., neglecting terms suppressed by X" , X', (p", (p'; hereafter prime 
denotes J-) are given in Appendix for completeness. However, we are interested in momenta 
P ~ e. According to the relations ([2]), modes of interest have relatively high spatial momenta 
and frequencies, 

P, > m, M , P, > H (4) 

where Q is the physical frequency. Note that in the regime ([31) one has H ~ M, so the 
second of these inequalities is a consequence of the first one. For momenta obeying (jlj), the 
equations simplify considerably: 

X" - A X + eab' - tadA - 2ea 2 X& = (5) 
eax' - 7 (&q - A&o) + 7aJf($" - A$) = (6) 
eadiX + a {b" - AbJ - 2aaXd^' = (7) 

-^a($" - A$) - 2eaX X ' + lX(b' Q ' - Ab ) + 2aXd i b' i = (8) 

47T 

We are now in a position to solve these equations in the WKB approximation by writing 



{x, b , bi, $} oc exp yi J cudr] — ipx 

with slowly varying amplitudes. Before doing that, we note that at P, Q ~ e, one estimates 
from eq. (IE|) 

Therefore, the last terms in eqs. ([5]), (EJ) and (JTj), describing back reaction of the gravita- 
tional potential on the field perturbations, are suppressed as compared to other terms by 



1 In fact, the behavior of modes in Minkowski space-time at P < M is somewhat more complicated |23) . 
This region of momenta is not of interest for our purposes, as it corresponds to physical momenta, and hence 
frequencies, small compared to the Hubble parameter, see 
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X 2 /Mp t ~ m 2 /e 2 . Neglecting these terms, we arrive at the system of equations for dark 
energy perturbations, which we write in the leading order of the WKB approximation, 

(uj 2 — p 2 )x — ieaubo + ieapbi = , 
ieaux + l(^ 2 - p 2 )b = , 
ieapx — a{uj 2 — p 2 )b L = 0, 

where b^ = [Pi/p)bi- This is basically the same system as that obtained in Minkowski space- 
time. There are three modes, two of which are normal at all momenta, while the third one 
is tachyonic at low momenta. In terms of physical frequency, the dispersion relation for the 
latter mode reads 

n 2 = ^ = p 2 + ^(i-Ji + ^P 2 \ . ( 9) 



where 

c = - + -. 

a 7 

As promised, the tachyonic regime occurs at 

p < Pc = 4= . 

The corresponding solution at P > P c is given by 

i(u 2 -p 2 ) i2L 
ecu 

bo = C-f(r})e^, 

7 

b L = - C .^f( v )e iK * , (10) 

where f{rj) is a slowly varying function, C is an overall constant and the phase factor is 
introduced for convenience. 

The function f(rj) can be found by noting that with our approximations, the action for 
the field perturbations reads 



drjdx. 



^m 2 - (d^) 2 ) + l{(b' Q ) 2 - (Vb ) 2 ) + U{ X ') 2 - (V X ) 2 ) + ea( X % - d lX h) 



(11) 

The corresponding energy functional is 



E= dx 



f ((^) 2 + (dibj) 2 ) + l((b' ) 2 + (V6 ) 2 ) + \((x r ) 2 + (Vx) 2 ) + ead lX h 
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In the WKB approximation, this energy is conserved. Making use of the expressions (TTUT) 
and requiring that energy of this solution is conserved, one finds 



(ui 2 —p 2 ) 2 
72 



(a 2 p 2 



(12) 



The same result can be obtained by considering the conservation of the Wronskian of the 
system (J5]), (jSJ), ©, which, with our approximations, is 



W 



-i (x*X + lK b 'o + ab l b> L + tax*b - c.c. 



Now, to find the overall constant in ([TO]) we quantize the system with the action (jlip at 
early times and obtain 

A + 

C = j2^> (13) 
where the creation and annihilation operators obey the standard commutational relation 

[^p,^] = <Kp-pO- 

We are interested in the behavior of perturbations in the tachyonic regime. The expo- 
nentially growing part is found by the standard WKB continuation of the expressions (jlOp . 
( I12p and (TT3|) "beyond the turning point". In this way we finally obtain in the tachyonic 
regime 



X 
h 

b; 



dp 



(2vr) 3 / 2 
dp 



i(uj 2 — p 2 ) 



/A+e" ipx + h.c. ) e 



JHdr] 



(2tt) 3 / 2 

where f(t]) is still given by (Tl2|) . 



(2l)3/2 ^/-4p + e- + h.c.je/M*, 



4 Gravitational potentials 

The gravitational potential $ is determined by eq. ([8]). After straightforward calculation we 
obtain that at P > P c the explicit expression is 

, /4tt m / 4 7 P 2 1 , 

*=-V275w — 7rW?TT 1 /( " )e 
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where C and f(rj) are the same as in (TlOT) . Continuing into the tachyonic region, we obtain 



An m 




dp 



27 M PL e J (2tt) 3 / 2 



4 7 p2 



1 + iC7!p2 _ ! 



- 1 /A+e- ipx + h.c. 



This expression gets simplified towards the end of the development of instability, i.e. at 
P <C e (but still P 3> H ~ M). In that case the dispersion relation ([9]) reads 



2 7 2 

CJ = p 

a 



so that 



Hence 



f(v) 
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3/4 a l/4 X 



(a + 7) 1 /2^ a ( ?7 ) • 



^_ t /g _jn 7 3/4 « 1/4 (a~7) /^P !_ M+p -*x 

V 27 Mpie (a + 7) 3 / 2 7 (2^)3/2^0 1 p 



+ h.c.) 



J |w|d»7 



This is Gaussian random field whose variance we write in terms of the physical momentum 
and frequency: 



where the moment t c corresponds to the beginning of the tachyonic regime, i.e., 

P 



(14) 



P{t c 



a(t c 



Note that the exponential growth factor is large but finite in the formal limit t — > oo. As an 
example, for constant H, the total growth factor is 



Q(t)dt 



1 

H 



dP e 
V(P)— = d(a, 1 )- 



H ' 



where d(a, 7) is of order one (for example, d(l, \) ~ 0.499). 



5 Discussion 



Recalling the relations (|2J) and ([3]) one observes that the gravitational potentials grow from 
very small values, set by vacuum fluctuations at the low momentum scale P c ~ e, but they 
eventually may become large. On the one hand, as discussed in Ref . [23J , this gives rise to a 
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constraint on the parameters of the model we discuss in this paper: by requiring that $ 1 
at the present time, one obtains 



where Hq is the present value of the Hubble parameter. So, the model may serve as a viable 
description of dark energy at the expense of extra fine tuning, over and beyond the usual 
fine tuning required to get the right value for cosmic acceleration. On the other hand, with 
even more fine tuning, the dark energy induced gravitational potentials may be sizeable (and 
still acceptably small) today. A novel feature here is that the emergence of these potentials 
would be a fairly recent phenomenon: the largest values of the exponential factor in (}T4"j) 
are obtained at late times when fl(t) changes slowly. Also, the amplitude of perturbations 
would be peaked at a certain momentum P: for high momenta, the tachyonic regime has not 
yet started, while for low momenta the integral in the exponent is saturated at early times, 
and hence is small. These properties may be potentially observable. 

The authors are indebted to O. Khovanskaya, M. Libanov and M. Sazhin for helpful 
discussions. This work is supported in part by Russian Foundation for Basic Research, grant 
08-02-00473. 

Appendix. Complete action and equations for pertur- 
bations 

For completeness, we present here the complete expressions for the quadratic action and 
linear equations for the scalar perturbations, valid in the slow roll regime. We use the 
conformal Newtonian gauge and neglect the terms suppressed by X", X', if", if'. The 
quadratic action is 




for a = 1 , 7 = 0.5 , 




(15) 
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where 



n V 8 2 2 8 V ; J 



a: 
2 



- aXdi® - aHbiY - a 2 H 2 b 2 + 2aH(dib - aXd^)b 



a 



- [(dibj) 2 - 2adibi(Hb - X&) + 3(-aHb + 4a'X$ + aX&) 2 ] 



-aa'X (di$bi - 3b & - 6a'X$ 2 - 6aX$$') + 



+ 
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(b' - aHb - aX&f - (<9;6 - aXd^f 



+ 



+- ( X ' - xa#) 2 - ~ (^X) 2 + ea (x' - X^H) (b - 2aX$) - ead iX h+ 



2 

+ — (6g - 4aX6 $ + 4a 2 X 2 $ 2 - o 2 ) - 



2^2 



ma 



-X • 



From here we obtain the equations for the field perturbations: 



X" - Ax X + m 2 a 2 x + eab' Q + 2ea'b Q - eadA - 2ea 2 X& - 6eaa'X$ = 

e{ax - a'x) - 7 (b^ - Ab - ^-b^j + (M 2 - 3aH 2 )a\ + 2aaHd i b l + 
+ 7 aX($" - A$) + 2(3a + 7)0! X& + 2(6a# 2 - M 2 )a 3 X$ = 

eada + 2aaHdib + a (b'! - Ah t - ^-b^j + (M 2 - aH 2 ) a\ - 2aaXd^' - 2aa 2 HXd^ = 

3 /If 2 

pl [ a 2^// _ + 2afl / $ / _ 2 ( aa // _ ( a ')2)$] _ 2eaX(ax' - a'x) + 7aX(^ - Ab )- 
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-6aa 2 HXb' - 7 a"X6 - 2M 2 a 3 Xb + 12aa 3 H 2 Xb - 6aa"Xb + 2aaXd i b' i + 
+a 2 X 2 (3a - 7)$" + (7 - a)a 2 X 2 A$ + 2(3a - 7 )a 3 M 2 $' + 6a(aa" - 5(a') 2 )X 2 $+ 

+4M 2 a 4 X 2 $ = 

Equations (EJ), ([6]), (J7J) and (JE]) follow from the latter equations in the regime (jlj), while the 
action (flBT) reduces to (ITT]) in this regime. 
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